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Abstract 

We analyze the asymptotic behavior of linear Fokker-Planck equations with 
time-dependent coefficients. Relaxation towards a Maxwellian distribution 
with time-dependent temperature is shown under explicitly computable con- 
ditions. We apply this result to the study of Brownian motion in granular 
gases as introduced in [7], by showing that the Homogenous Cooling State 
attracts any solution at an algebraic rate. 
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1 Introduction 

The linear Fokker-Planck equation (FPE) 

^(v, t)=XV- {v(p(v, t) + 6Vcp(v, t)} , v E R N {N ^ 1) (1.1) 

arises in many fields of applied sciences such that statistical mechanics, chemistry, math- 
ematical finance (see the monographs |19j and ^U] for a large account of applications to 
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the FPE). Such a drift-diffusion equation can be derived from the Langevin equation to 
model the Brownian motion of particles in thermodynamical equilibrium. In this case, the 
parameters A and 9 are two positive constants which represent respectively the friction 
term and the temperature of the system. The qualitative analysis of equation 1)1. 1|) is well 
documented in the literature. We refer to ^Hl f° r a precise description of the hilbertian 
and spectral methods used in the study of 1)1.1(1 and to [21 IE] for the relatively recent 
approach to the L 1 -theory by means of entropy-dissipation methods. 
It is easy to notice that the set 

M = J g G L 1 ^), g ^ 0, / g(v)dv = l, / vg{v)dv = 0, / v 2 g(v)dv = 9 < oo I 
I Jm Jr n Jr n J 

is invariant under the action of the right-hand-side of p. 1(1 . Moreover, it is well-known that 
(|1.1|) admits a unique steady state in M given by the Gaussian distribution (Maxwellian 
function in the language of kinetic theory) 

M e {v) = (2^)- Af / 2 exp {-v 2 /29), v G R N . 

Entropy-dissipation methods (see [3] for a review on recent results on the topic) provide 
a precise picture of the asymptotic behavior of the solution to (|1.1[) for initial data in M. 
Given / £ M, the (Boltzmann) relative entropy (finite or not) of / is defined as 

H(f | M f ) = H(f) - H(Mf) = [ /log (Jfr) dv 

Jm. n \ M fJ 

where Mf is the unique Maxwellian distribution in M with the same temperature as /. 
Given /o 6 M, with the assumption of bounded initial relative entropy 

H(f o \M ) < c>o, 

it has been proven in [2j that the unique mass-preserving solution f(v,t) of (|l.lj) decays 
exponentially fast with rate 2A to Mq(v) in relative entropy, i. e, the estimate 

H(f\M e )(t) <e- 2Xt H(f \M e ) 

holds. The classical Csiszar-Kullback inequality then allows to translate convergence in 
relative entropy to the more standard L 1 -setting. The following result is proved in [2]: 

Theorem 1.1. Let A and 9 be two positive constants. Let us assume that /o £ M has a 
finite relative entropy. Then, there exists a constant C > depending only on the initial 
relative entropy such that the solution f(v,t) to (|1.1|) fulfills 

\\f(;t)-M e \\ L i {RN) ^Cexp(-\t) 

for any t ^ 0. 
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Main objective of this paper is to generalize this result allowing the friction term A 
and the temperature to fluctuate with time. In this case, the Fokker-Planck equation 
reads: 



where X(t) and 6(t) are positive functions of time. 

Non-autonomous Fokker-Planck equations arise for instance in the study of a periodi- 
cally driven Brownian rotor and in this case \(t) and 9{t) are periodic functions of time. 
In statistical mechanics, equation (|1.2|) arises as a natural generalization of equation 
in the context of non- equilibrium thermodynamics Among other models, equation 

(|1.2|) appears in the study of the tagged particle dynamics of a heavy particle in a gas of 
much lighter inelastic particles. As observed by J. J. Brey, W. Dufty and A. Santos 0, 
the large particles exhibit Brownian motion and the Boltzmann-Lorentz kinetic equation 
satisfied by the distribution function of large particles can be reduced to a Fokker-Planck 
equation whose coefficients depend on the temperature of the surrounding gas. Granular 
gases being non-equilibrium systems, this temperature turns out to be time-dependent 
and the Fokker-Planck equation derived in [Jj is of the shape (|1.2|) . Since the study of the 
long-time behavior of the solution to the Brey-Dufty-Santos model is one of the main 
goals of our analysis, we will explain with much more details the approach of [Zj in the 
next section. 

Because of the time-dependence of both A(i) and 6(t), equation (|1.2|) does not possess 
stationary states. Nevertheless, two natural questions arise: 

- Do they exist particular (self-similar) solutions to (|l,2j) which attract all other so- 
lutions (as the Maxwellian does in the autonomous case)? 

- If such self-similar solutions exist, is it possible to reckon the rate at which they 
attract the other solutions? 

We answer positively to these two questions under some reasonable conditions on the 
time-behavior of friction and temperature. Our method is based upon suitable (time- 
dependent) scalings which allow us to transform the non-autonomous equation 1)1.2(1 into 
a Fokker-Planck equation of the form (jl.l|) and then to make use of Theorem 11.11 Clearly, 
the self-similar profile is a Maxwellian with time-dependent temperature. In this context, 
the Maxwellian distribution plays the role of the Barenblatt profile in the study of the 
porous medium equation [2T] . 

Application of this abstract result to our motivating example in kinetic theory of 
granular gases shows that the distribution of the Brownian particles relaxes towards a 
Maxwellian distribution with time-dependent temperature. While this fact has already 
been noticed by J. Javier Brey et al. [Zj our result gives a precise estimate of the rate of 
convergence towards this self-similar profile (the so-called Homogeneous Cooling State) 
which turns out to be only algebraic in time. Analogous results, which try to clarify the 
role of the Homogeneous Cooling State in kinetic models of granular gases, have been 
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recently obtained for the case of the Boltzmann equation for inelastic Maxwell particles 
jl]. We postpone a detailed discussion on this point in the conclusions of this note. 

The organization of the paper is the following. In the next section, we present in 
some details the derivation of the FPE in the context of Brownian motion for granular 
gases. In Section 3, we deals with a general non-autonomous FPE and we answer to the 
two aforementioned questions (Theorem I3.1j) . Finally, in Section 4 we turn back to our 
motivated example of Brownian particles and we show how the abstract result of Section 
3 allows us to estimate the rate of convergence towards the Homogeneous cooling state. 



2 The Brownian motion in granular gases 

The motion of heavy granular particles of mass m embedded in a low density gas 
whose particles have mass m g with m g <C m has been considered in [7]. The particles 
under consideration are assumed to be hard-spheres of IR 3 and, for the sake of simplicity, 
the diameters of the particles of both species are assumed to be equal and normalized 
to unit. The case of particles with different diameter can be investigated as well, and 
does not lead to major supplementary difficulties [Jj. The collisions between the heavy 
particles and the fluid ones are partially inelastic and are characterized by a coefficient 
of restitution e £ (0, 1). Assuming that the concentration of heavy particles is small, one 
neglects the collision phenomena between them. Let us denote by f(v,t) the distribution 
function of the heavy particles having velocity v £ M 3 at time t > and by g(v,t) the 
distribution function of the surrounding gas where, for simplicity, it is assumed that these 
two quantities are independent of the position. Then, the evolution of f(-,t) is given by 
the Boltzmann-Lorentz equation, which in weak form reads 

d7 



f(v,t)i/>(v)dv= \q-n\f(v,t)g(w,t)ty(v*)-ip(y)]dMwdn (2.1) 

'R 3 xR 3 xS 2 



for any test-function ip(v). Here q = v — w and (v*, w*) are the post-collisional velocities: 

A(l + e), A(l + e), , x 

v*=v- ^ ' (q-n)n, ^ = «, + _A_Z( q . n )n (2.2) 

where A is the mass ratio A = m g /m. Note that, by assumption, A <C 1. 

To solve the linear equation (|2.1|) . one has to make explicit g(v,t). Assuming that 
the binary collisions between the surrounding particles are inelastic and characterized 
by a constant restitution coefficient < e g < 1, g(v,t) is given by a solution to the 
(nonlinear) Boltzmann equation for granular hard-spheres jSl EH- Leaving details to 
the pertinent literature, the role of the equilibrium Maxwellian function in the elastic 
Boltzmann equation is here represented by the Homogeneous Cooling State (see References 
[TT| IT3] ) which implies that 

g(v,t) = v g (t)- 3 $ 



Vg{t) 
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where v g (t) is the thermal velocity of the gas particles defined as v g (t) = [2T g (t)/m g ] ' . 
The temperature T g (t) is defined in a standard way (see Section 3). The self-similar profile 
<£(•) is a stationary solution of some suitable steady Boltzmann equation (see ^DIH]) and 
is not explicitly known. However, an important fact to be noticed is that &(v) is a 
function of e g which, in the quasi-elastic regime e g — > 1, converges toward the Maxwellian 
distribution 7r -3 / 2 exp (— v 2 ). The temperature T g (t) is cooling because of the inelasticity 
of the collisions. Hereafter, we will assume T g (t) to obey the the so-called Haff's law [Tlj : 

T g (t)=T g (0)(l + t/r )- 2 (2.3) 
where To > is the characteristic time 



T n 1 = — — e A f <$>(v)<S>(ti)\i-w\ 3 dvdw. (2.4) 

12 y 2m g J R s xR 3 

The Boltzmann-Lorentz equation ()2.1|) can be reduced to a FPE with time-dependent 
coefficients of the form (|1.2j) performing two asymptotic procedures: 

• The first procedure is a simple extension of the standard method for elastic particles 
(grazing collisions asymptotics [U], see also |16| in the context of the dissipative linear 
Boltzmann equation). Precisely, according to (|2.2|) . one sees that, when a heavy particle 
collides with a small one, the velocity of the heavy particle is only slightly altered: 



A(l + e), 

'q • n)n 



1 + A 



< 1 



so that v* ~ v. Therefore, performing in 1)2. 1|) a formal expansion to leading order in 
the mass ratio as A — > 0, one obtains the following Fokker-Planck equation with time- 
dependent coefficients: 



df(v,t) _ 

- Vjj 



A(v,t)f(i,t) + ■ (N(v,t)f(v,t)) 



dt 

where the vector A(u, t) and, respectively, the tensor N(y, t) are given by 



(2.5) 



A ._ . (1 + e Avr f . , (J „ 



and 

N ^ 5, ^ = (l+r) ^^ 3 9(*>*) (|q| 3 ^' + 3|q|qiqi)du; (i,j = 1,2,3). 
We used above the standard notation 




V,-(N(v)f(v))= \ Y—N i:j (v,t)f(v,t) 



i=l,2,3 
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C = |^(l + e)A / $(w)\w\dw, (2.7) 



(see |7J Appendix A] for a detailed derivation). 

• To simplify further the Fokker-Planck equation 1)2.5 J) one performs a second asymp- 
totic procedure which consists in assuming that the thermal velocity of the heavy particles 

1 /2 

v(t) = [2T(t)/m] 1 is negligible with respect to the one of the surrounding particles v g (t): 
v(t) <C v g (t). This leads to a formal expansion in T(t)A/T g (t) where T(t) is the tempera- 
ture of the Brownian particles (see Remark l2.1)) . In this case, the vector A(v,t) and the 
tensor N(v, t) reduce to 

A(v,t) ~ a(t)v and Nij(v,t) ~ 2rj(t)Sij, i,j = 1,2,3 

where 

a(t) = (Tgii) 1 ? 2 and rj(t) = iT g {tf' 2 . (2.6) 
In 1)2.6)) we used the notations 

2\/2vr 
and 

-3/2 . 

C = ^V(l + e ) 2 A 2 / Hw)\wfdw. 
3V2 Jr3 

Taking these simplifications into account, equation Q2.5JI is replaced by the following 
df 

-±(v, t) = V e ■ {a(t)v f(v, t) + V (t)Vef{v, t)} . (2.8) 

The Fokker-Planck equation 1)2.8)1 is of the form ()1.2)1 . The drift and diffusion coefficients 
depend on time only through the surrounding gas temperature T g (t) (see 1)2.6)1 ). 

Remark 2.1. The derivation of ()2.8)) from 1)2.1)1 is based upon a time-dependent asymp- 
totic procedure where it is assumed that 

It- 

As shown in [2j this assumption requires A —* and e 9 — > 1 (quasi-elastic regime for the 
surrounding gas) simultaneously. A further consequence of these assumptions is that 

9 < 1. (2.10) 



2 V / 2A 1 

We will find again this condition hereafter. 
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3 Long time behavior of non- autonomous Fokker— 
Planck equations 

In this section we consider the general Fokker-Planck collision operator with time- 
dependent coefficients written in the divergence form: 

Q t (f)(v)=\(t)V v -{vf(v,t) + 8(t)V v f(v,t)}, v£R N (N^l) (3.1) 

where X(t) and 9{t) are two positive functions of time. We are concerned with the large- 
time asymptotic behavior of the solution to the Cauchy problem 

^(v,t) = Q t (f)(v) veR N ,t>0 (32) 
f(v,0)=f (v) 

where the initial data /o is assumed to be nonnegative and integrable, 

f >0, and / € L 1 (R JV ). 

In accordance with the language of kinetic theory, we define the mass density g(t), mean 
velocity u(i) and temperature T(t) respectively as: 

£>(*)=/ f(v,t)dv, u(t) = -7-r[ vf(v,t)dv, 

and 



Ng{t) J r n 

The number density is preserved by the (non-autonomous) Fokker-Planck operator 
while the mean velocity is preserved only if initially equal to zero. Precisely, if f^ N vfo(v)dv - 
0, then u(t) = for any t > 0. In this case, the evolution of the temperature is 

= _2A(t) (T(t) - 6(t)) , (t>0). (3.3) 

In order to find the intermediate asymptotic for ()3.2j) . we look for a solution to (|3.1|) of 
the shape: 

f(v, t) = j(t)- N F (vh{t), r{t)) = j(t)- N F(v, t) 

where the new time scale r = r(t) is nonnegative and such that r(0) = 0, the scaled 
velocity is 

v = v/^(t) 

and 7(-) is positive. Without loss of generality, one may assume that 7 2 (0) = To := 
Jr3 v 2 fo(v)dv so that 

F(v,0) = F(v/T ,0) = fo(v/T ). 



7 



One sees immediately that 



df f(t) dF 



(3.4) 



where the dot symbol stands for the time derivative. In the same way, one can show that 

(3.5) 



= JfV, • (vF(v,r)) + ^lvlF(r. T ). 



This leads to the following evolution equation for F(-,r): 

7(*) 



1 

^' T) = ^) 



A(t) + 



7(*) 



V s - (vF(v,t)) + 



X(t)8(t) 2 
rWW ViF(v) 



One notes that (|3,6|) reduces to a "good" Fokker-Planck equation 
<9F 



(5, r) = Vfi • (C r) + Vei^w, r)) u G 1^, r > 



dr 



provided 



and 



f(t) 



A(i) + 



X(t)0(t) 



i 



1 Vt > 0, 



Vt > 0. 



(3.6) 

(3.7) 

(3.8) 
(3.9) 



r(t)7(*) 2 

Of course, to investigate the asymptotic behavior of F(-,t) and apply Theorem 11.11 one 
has to find conditions on A(-) and #(•) insuring that the time scale r verifies 

lim r(t) = +oo. 

t^oo 

Solving equation (|3.8|) - (j3.9|) leads to 

X(t)9(t) 



I 2 it) 



f(t) = A(i) + 



m 

7(0 



i. e. 



1 d 



X(t)0(t) = 7 2 (t)A(i) + 7 (t)7(t) = l\t)X(t) + - ^{7 2 WI- 
Since 7(0) = \/7q, one obtains 

7 (t)=exp^-^ A(s)ds^) |r + 2^ A(s)#(s) exp ^2 ^ X(r)dr^j ds 



1/2 



t > 0. 
(3.10) 



S 



Now, from (f3~3j). 



f(t) 



1 d i 



\(t)6(t) exp ^2^ A(s)ds 



T + 2 A(s)0(s) exp ^2 A(r)dr^ ds 
T o + 2 J q A(s)#(s)exp ^2^ A(r)dr^ds 



Solving this equation with the initial datum r(0) =0 one gets 



r(t) = \ log 



T + 2^ A(s)0(s)exp A(r)dr^ ds 



Clearly, 



lim r(i) = +oo if and only if / A(t)0(i) exp I 2 / 



(i>0). 



A(s)ds I dt = oo. 



(3.11) 



This leads to the following result. 

Theorem 3.1. Let us assume that A(-) and9(-) are nonnegative functions onM + satisfying 

r-t 



j \{t)9(t) exp ^2 j A(s)ds^) dt = oo. 



(3.12) 



Let us assume furthermore that /o E M has a finite relative entropy. Then, there exists a 
constant C > such that 

\\f(;t)-foo(;t)\\Li(W) < ^ *>°- 



T + 2 A(s)#(s) exp 2 / A(r)dr ds 



1/2 



(3.13) 

The intermediate asymptotic profile foo(v,t) is given by 

/«,(«, t) = (27TT(t)y N / 2 exp{-v 2 /2T(t)} = M T(t) (v) 
where T(t) is the temperature of f(-,t) given by (j3.1U|) , 

Proof : The proof reduces to the study of (|3,6|) . Clearly, one may choose a = 1 (this is 
equivalent to change r by (|3.11jl ). By (|3.12|) . r(t) — > oo, and according to Theorem ll.il 

\\F(.,T)-M 1 (.)\\ Ll{R N ) ^Ce W {-T} r>0 (3.14) 

provided F(v, 0) is of finite relative entropy. This is the case since F(v,0) = fo(v/To). 
Now, turning back to the original variables, one gets the conclusion using the fact that 



exp{-r(t)} = |l + 2^ A(s)(9(s)exp ^2^ A(r)drJ ds 
by virtue of p.llj) . 



-1/2 



(t>0) 
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Remark 3.2. The intermediate asymptotic is given by the Maxwellian distribution with 
the same temperature T(t) as the one of /(-,£)• Of course, this Maxwellian distribution 
is a particular solution to (|3.2|) . The most important feature of Theorem 13.11 is that it 
provides the rate of convergence of any solution f(v,t) towards the self-similar profile 
foo{v,t). This rate is explicit in terms of the known coefficients X(t) and 0(t). 



4 The homogeneous cooling state for the Brown- 
ian particles 

We apply here the results of Section 3 to the study of the so-called homogeneous 
cooling state for equation (|2.8|) . Let fo(v) be an element of M and let us consider the 
Cauchy problem 

^(v,t) = \(t)Vv-{vf(v,t) + e(t)Vzf(v,t)} veR 3 ,t>0 
f(v,0) = f (v) 

where we transformed the right-hand side of equation (|2.8|) into a non-autonomous Fokker- 
Planck operator of the form (|3.1|) by setting, according to (|2.6|) . 

A(t) = a(t) = QT g {t) 1 / 2 and 9{t) = 4| = iC%{t). (4.2) 

Accordingly, for any t > and with the notations of Section 3 

J \(s)d{s)exp (2 \{r)dr\ds = £ j T g (sf /2 exp ^2C J ^T g (r)dr^J ds 

= CT g (0) 3/2 J\t + s/r )- 3 e W fe^O) J' T^r) ds - 
Note that the last equality follows from Haff's law ()2.3|) . Clearly, for any s > 

(i+ s /r r 3 exp fey/rM^T^^) = a+^or 



with v = 2(roy / T g (0) — 3. Consequently, condition (|3. 12|) of Theorem 13.11 is verified pro- 
vided v ^ — 1. Note however that, if v = — 1, then 

T(t)=exp^-2^ \{s)ds^ |t + 2^ X(s)9(s) exp ^2 J \{r)dr^j ds\ 
= (1 + t/ro)' 2 (Tq + £T 9 (0) 3 / 2 t log(l + t/r )) , 

so that 

T(t) 

— . . — > 00 as t — > 00. 
10 



In this case, assumption (|2.9|) is violated, and the Fokker-Planck equation (|2.8[) is mean- 
ingless. Hence, condition ([3.12ft of Theorem 13 .11 reduces to v > —1, i. e. 



(r ^T g (0) > 1. (4.3) 
Using (|2.4|) and (|2.7|) the above condition reads 

i6 (l + 6)A f$(v)\v\dv > 1 

l — €g f$(v)$(w)\v - w\ 3 dvdw 

Now, recall that the Fokker-Planck equation ([3.1[1 turns out to be valid only for nearly 
elastic surrounding particles (see Remark 12.1(1 . Since in this quasi-elastic regime <&(w) 
approaches the Maxwellian distribution tt~ 3 ^ 2 exp(— v 2 ), one can reasonably approximate 
the two moments f R3 $>(v)\v\dv and f R 3 x ^3 $(v)$(w)\v— w\ 3 dvdw by their limits as e g — ► 1, 
obtaining 

/ 7r _3/2 exp (-v 2 )\v\dv = 2/yJ~n 
it 3 

and 

/ 7r~ 3 exp (— v 2 — w 2 ) \v — w\ 3 dvdw = \Q/V2t:. 
Using these approximation, equation (|4,3|) turns out to be equivalent to 

1 — e 

9 

Once again, we find the same condition (|2.1UI) of validity of the non-autonomous Fokker- 
Planck equation ([2.8)1 . 

Remark 4.1. Note that, here again, we only assume that the surrounding gas particles 
suffer nearly elastic collisions (i.e. 1 — e g <C 1) but we not assume e g to be equal to 
one. As a consequence, we do not replace the cooling state profile $ by the Maxwellian 
distribution, but we only assume that its moments do not differ to much from the ones of 
the Maxwellian distribution. 

The previous reasoning leads to the following Theorem. 

Theorem 4.2. Let us assume that 

— = 9 - — ► /T 1 < 1 as A 0, t Q -> 1. 

2 V / 2A 1 + e 9 

Let /o G M be of finite relative entropy. Then, the solution f(v,t) to the Fokker-Planck 
equation ()3.1[) converges towards the cooling Maxwellian 



/«,(«,*) = (2vrT(t))- 3 / 2 exp{- ? ) 2 /2T(t)} , 
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and the following bound holds 

\\f(;t) - /oo(-,t)||n(R3) = O(^) (* - oo). 
The temperature of the Maxwellian is given by 

_ (l + 6)AT g (0) 2 / (l + e)AT 9 (0) \ 2/3 

r(t) - 2m fl (i - (1+ * /ro) H /o) ' 

Proof : The proof is a straightforward application of Theorem IH. II Here we use the fact 
that 

J A(s)0(s)exp (2 J A(r)dr^ds = ^T 9 (0) 3/2 ^ (1 + s/r ) 2/3 " 3 ds 

= 2^2^ (0)3/2 { (1 + ' /T0)2 ^ 2 - 1 }- 



Moreover, by (|3.10f) . the temperature of the heavy particles is 

T(i)=exp^-2^ A(s)ds^) |r(0) + 2^ A(s)0(s) exp ^2 jf A(r)drjdsJ. 
This last quantity is explicitly computable using the expressions of tq and £. ■ 

Remark 4.3. One notes that T(t) obeys asymptotically Haff's law since the decay of 
temperature of the heavy particles for large t is in O ((1 + t/T )- 2 ). Actually, 



T g (t) 2m,(l -/J- 1 ) 

An interesting feature is that, depending on the values of the parameters e, A, m g and 
f3, the temperature of the heavy particles is greater or smaller than the one of the sur- 
rounding gas. This contrasts the classical case of elastic particles in equilibrium. Indeed, 
in this case, according to Theorem ll.il the distribution function relaxes to a Maxwellian 
distribution whose temperature 6 is exactly the one of the surrounding bath (see [2] for 
more details). We refer the reader to [7j for a discussion of the competing effects which 
imply the asymptotic difference between T(t) and T g (t). 

Remark 4.4. We point out that, whereas for systems in equilibrium, the relaxation rate 
is exponential (Theorem I l.ljl . one notes here that the Homogeneous Cooling State foo(v,t) 
attracts the distribution function f(v,t) only with an algebraic rate. 
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5 Concluding remarks 



We discussed in this paper the intermediate asymptotics of a linear Fokker-Planck 
equation with time-dependent coefficients of the form Q1.2|) . We showed that, under some 
reasonable assumptions on the drift and diffusion coefficients, any solution f(v,t) to (|1.2jl 
relaxes towards a Maxwellian distribution function whose (time-dependent) temperature 
is the one of f(v,t). More important is the fact that the rate of convergence towards this 
self-similar solution is explicitly computable in terms of the coefficients and the initial 
temperature. 

We applied our result to the motivating example of Brownian motion in granular 
fluids, already addressed in [7]. For such a model, the Fokker-Planck equation (|1.2|) is 
an approximation of the Boltzmann-Lorentz equation. According to our general result 
( Theorem 13. the so-called Homogeneous Cooling State for this model is a Maxwellian 
distribution whose temperature obeys asymptotically the Haff 's law. Moreover, the rate of 
convergence towards this self-similar solution is algebraic in time. We wish to emphasize 
here the fact that the question of the rate of convergence towards the Homogenous Cooling 
State (HCS) is of primary importance in the kinetic theory of gases. We recall here that any 
solution to the non-linear Boltzmann equation for inelastic interactions relaxes towards a 
Dirac mass because of the dissipation of the kinetic energy. It has been conjectured however 
by Ernst and Brito [111 I12j that the HCS attracts any solution faster than the Dirac mass 
does. While for hard-spheres interactions, only few results support this conjecture, the 
situation is almost clear for the non realistic case of the Boltzmann equation for maxwellian 
molecules jl]. For this model in fact, it has been shown that the (unique) self-similar 
solution attracts any other solution corresponding to an initial density with a sufficiently 
high number of bounded moments. Other simplified models have been studied recently. 
For nearly elastic flows in one-dimension, the nonlinear Boltzmann equation reduces to a 
nonlinear friction equation |18[I20| and it has been shown in (Sj that, in this case, the HCS 
does not attract much faster than the Dirac mass since the improvement in the rate of 
convergence is only logarithmic in time. This question has also been addressed recently in 
|15j for general nonlinear friction equations corresponding to a relative velocity dependent 
coefficient of restitution. 
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